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We consider a model of N two-dimensional bosons in a harmonic potential with weak repul- 
sive delta-function interactions. We show analytically that, for angular momentum L < N, the 
elementary symmetric polynomials of particle coordinates measured from the center of mass are 
exact eigenstates with energy N(N — L/2 — 1). Extensive numerical analysis confirms that these 
are actually the ground states, but we are currently unable to prove this analytically. The special 
case L = N can be thought of as the generalisation of the usual superfluid one-vortex state to 
Bose-Einstein condensates in a trap. 



I. INTRODUCTION 

The recent observation of Bose-Einstein condensa- 
tion in dilute gases of alkali metals JijH has stimulated 
much interest in the properties of systems of interacting 
bosons. A question of particular interest is whether such 
systems will form vortices under rotation, as occurs in 
superfluid 4 Hc. Experimentally such vortex states have 
been observed, both in two component systems B and 
in a stirred condensate of 87 Rb atoms || . Theoretically 
the stability of such vortices has been considered both in 
the Thomas-Fermi limit of strong interactions between 
atoms using the Gross-Pitaevskii equation HQ , and also 
in the weak interaction limit [p[-p5[ where the coherence 
length is much larger than the size of the atom cloud. It 
is the latter limit we shall focus upon in this paper. 

Pursuing an analogy with the fractional quantum 
Hall effect B, we previously introduced a model of 
weakly interacting bosons in a harmonic well || to ad- 
dress the question of whether attractive bosons condense. 
We proved analytically that all the angular momentum 
in this model is carried by the center of mass for attrac- 
tive bosons, whereas for repulsive bosons we numerically 
found that a vortex state forms in the case of one unit of 
angular momentum per boson. Further numerical work 
by one of the authors [ p~0|JTl] ] extensively investigated the 
properties of ground states of the repulsive model for an- 
gular momentum L > N , demonstrating that although 
such states are more complicated than the analytic ones 
known for L < N, they can still be understood either 
within vortex or composite fermion or boson pictures. 
Mottelson |l2| considered the low-lying eigenstates for 
the case L < N, whilst Bertsch and Papenbrock |l5| ] 
performed numerical computations and noted that the 
ground state of the repulsive model for L < N is the 
elementary symmetric polynomial, §£,, of coordinates, 
zi = Xi + iyi, relative to the center of mass, R — Zi/N . 
Finally, recent work by Kavoulakis et al [jlij and Jackson 
et al |lq | has considered the relationship between mean 
and exact numerical solutions in the limit of large N . 

In this Letter we present an analytical proof that 
the state described above is an exact eigenstate of the 
model for L < N. We have unfortunately been unable 
to show that this state is the ground state, although we 



have considerable numerical evidence to show that this 
is the case. 

The model is of N bosons in a harmonic potential 
in two dimensions interacting via a delta-function poten- 
tial, for which the Hamiltonian is 
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The natural way to look at this problem is in second 
quantised form, so first we must solve the non-interacting 
problem. We do this in plane polar coordinates since 
we are interested in angular momentum properties. The 
one-particle wavefunctions and energies arc then 
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where R(n, l,x) is the confluent hypergeomctric function, 
n r is the radial quantum number and t is the angular mo- 
mentum. The total energy and angular momentum of a 
system of N non-interacting bosons in this harmonic well 
is thus 
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The ground state manifold is then obtained by putting all 
bosons into the lowest radial state, n r i = 0, and choos- 
ing angular momentum ti all of the same sign (which we 
choose to be positive), such that 



Y,^ = L E tot = {L + N)hLo. 



(4) 



It can be seen that the ground state manifold has a de- 
generacy which is equal to the number of partitions of the 
total angular momentum L into the TV angular momenta 
ti of the individual bosons. 

If we assume that the dimensionless interaction 
strength is very small, |ry| <C 1, where rj = V/fiu), we 
can treat it as a perturbation whose sole effect will be 
to lift the large degeneracy of the non-interacting ground 
state. This means we can work within the non-interacting 
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ground state manifold, which freezes out the one-particle 
part of the Hamiltonian, and consider the effect of the 
interaction using degenerate perturbation theory. This 
is essentially an extension of an approach used for frac- 
tional quantum Hall effect systems Q to bosons. The 
normalised one-particle wavefunctions are given by 
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where we have moved to complex notation. The second 
quantized form of the interaction Hamiltonian is thus 



H ^ ' V mim2nin2 a mi a^ n2 a ni a n2 



mi ,m2,ni ,712 

where the matrix element 



III i //<•_>// ; 



r) (mi + n±)\ 
47T 2 mi+m2 



J mi+m2,rii+n2* 



(0) 
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For future convenience we will set r/ = 4ir. 

To perform the degenerate perturbation theory for 
given N and L, we first need the basis states which are 
labelled by the partitions of L into N pieces. Let us write 
a partition A in the form 



A = o Ao l Al 2 A2 
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where ^2 Xi — N and ^iXi = L. The corresponding 
basis state is then 
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where creates a boson of angular momentum I. In 
coordinate space, this basis state takes the form 
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where m\ is the monomial symmetric polynomial corre- 
sponding to the partition A. The latter is the symmetric 
polynomial of the N variables (z\, z% . . . z/v) which has Xi 
i-th powers. 

To consider the problem either analytically or nu- 
merically requires us to calculate the elements of the 
symmetric interaction matrix, — (X\V\fi). This is 
obviously best performed using the second quantized ap- 
proach. There are two types of matrix elements, diagonal 
and off-diagonal, and their evaluation is different. For 
the diagonal elements H\\ , we must sum over every pos- 
sible pair of elements in the partition, whether distinct 
or identical, 
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The off-diagonal elements can only be non-zero if A 
and n differ only in the angular momenta of two parti- 
cles. This can affect either 4 separate angular momen- 
tum states, as in the case where the angular momentum 
transfer is + 4— >l + 3;or3 separate angular momen- 
tum states, as in the case of angular momentum transfer 
+ 4^2 + 2. In terms of partitions, in case (i) Xi = fa 
except at four values of i, which we call and 
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The matrix element is then 



— H^x — A^/Xi 1 Xi 2 fii 3 /Xj 4 Vi x i 



(12) 



In case (ii) A^ = fii except at three values of i which we 
caUii,i 2 ,« 3 , and X il = X i2 =/x <2 -l, A, 3 =/i ia -l. 

The matrix element is then 



H\n — — 2y/Xi 1 (Xi 1 — 1) fa 2 fa 3 Vi 1 i 1 i 1 



(13) 



These formulas allow one to write down the interaction 
matrix Hxfj, which can then be diagonalised to give the 
energy eigenvalues and eigenstates. 



II. THE SUBSPACE PROPERTY 

If we look at the original Hamiltonian, we find that 
the interaction term depends only upon relative coordi- 
nates. To see this, change variables to the center of mass 
variable, R = J2i z i/N, and N — 1 relative coordinates 
such as Zi = Zi — R, where i = 1 ... N — 1: the inter- 
action is then a function only of the relative coordinates 
Zi. It follows that if ip{ z i ■ ■ - z n) is an eigenfunction of 
H with a certain energy E and angular momentum L, 
then Rij)(z\ . . . zn) is an eigenfunction of H with energy 
E and angular momentum L + 1. This means that all 
eigenf unctions of H at a given L are reproduced at all 
higher L, and thus the number of new states at any L 
is just til — riL-i, the difference between the number of 
partitions of L and L — 1 respectively || . 

This subspace property makes it natural to think in 
terms of a second type of symmetric polynomial, the el- 
ementary symmetric polynomials. These are defined by, 
for N variables, 



e L = ^2 z n z l2 ...z iL (14) 

il<i2<---<iL 

where L < N. For a general partition A we define 



ex = lie 



Xi 
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The set of new eigenstates at any total angular momen- 
tum L is seen to be spanned by the ex, the elemen- 
tary symmetric polynomials of the relative coordinates 
Zi, where we now include zn — —Z\ — Ufa... Zff-i- Since 
ei = 0, only partitions with Ai = can be formed (i.e. 
partitions with no part equal to 1), which gives exactly 
the correct number of states. 
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III. PROOF THAT THE E L ARE EXACT 
EIGENSTATES 



The final result for the operation of the Hamiltonian on 
ex, is thus 



In this section we prove that the states ex for 
L < JV (L / 1) are eigenstates of H with eigenvalue 
JV(JV — 1 — L/2). The proof is essentially a brute force 
method: we operate the Hamiltonian on the state ex, 
and show that the result is the above eigenvalue times 
§£. The derivation naturally falls into five steps, which 
we detail below. 

(1) Writing ex in Terms of eL and R 

Consider e^/, which can be written as 



R). 



(16) 



*1<*2<.--<«A/ 



If we expand out this product we will get the ele- 
mentary symmetric polynomials of the z,, namely the 
ex,, L < M, multiplied by R M ~ L . To get the cor- 
rect coefficients in this expansion we note that §m has 
AT! /Ml (JV - M)\ terms whilst e L has JV! /LI (JV - L)l 
terms. In the expansion of eu, each product term will 
produce Ml /LI [M—L)\ terms of the type which will add 
up to produce R M ~ L eL, so that the coefficient of ex in 
the expansion is (JV - L)\ /(JV - M)l (M - L)\ It follows 
that 



M 
X=2 



M — L N - LI jjM-L 



JV - M ! M - L! 



+( _ 1)A /-i^!(^-l) D Af 



JV - Ml Ml 



where we have noticed that, since e\ = R, the last two 
terms have the same form and should be combined. 

(2) Operating Hamiltonian H on eL 

An important feature of ex is that it is also the 
monomial function m\ corresponding to A = N ~ L 1 L . 
The normalised version of ex can thus be written as 
\e L ) = \0 N ~ L l L ), and it is clear that the Hamilto- 
nian H can only connect this state to itself and | A ) = 
| JV ~ i+1 l L-2 2 1 ), where we have labelled this state as 
| A ) for convenience in what follows. The two matrix el- 
ements can be calculated using from the formulas derived 
in the previous sections for diagonal and off-diagonal el- 
ements respectively. The diagonal element is given by 



(e L \H\e L ) ==N 2 -N- -L(L - 1). 



(18) 



The off-diagonal element is found by using the rule for 
the case where only three separate angular momentum 
states change (here 1 + 1 — > + 2) to give 
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■A\H\e L ) = -y/2L(L-l)(N-L- 
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H\e L ) = 



N 2 - N - -L(L — 1) 



\e L ) 



+ - y/2L(L - 1)(N - L + 1)\ A) . (20) 

(3) Removing the Normalisation Factors 

We want to get rid of the normalisation factors for the 
eigenstates so that the Hamiltonian will act directly on 
symmetric polynomials such as ex . In the previous equa- 
tion we should divide by the normalization factor for 
\0 N - L 1 L ) and multiply by that for | Ar ~ i+1 l L ~ 2 2 1 ). Us- 
ing the form for the normalisation factors given in Eq. 
( |io| ) , we find that the normalised version of Eq. ( po| ) is 

H\e L ) — [JV 2 — JV — \L[L - l)\\e L ) + — £±1 \A), 

(21) 

where the |A) are the symmetric polynomials with no 
normalisation factor. 

(4) Relating | N - L+1 1 L - 2 2 1 ) to e L and Re L _i 

Consider the product 



NR\e L -i) = 



' N 

_i=l J ii<. ..<■>£_! 
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(17) The above product can clearly only produce \0 N ~ L 1 L ) 
or |0 Ar_L+1 l i_2 2 1 ), depending upon whether the Zj 
from the prefactor is not or is included in the set 
(zi i: Zi 2 . . . Zi L _ 1 ). Moreover we see that each element of 
|0 W_L+1 1 L_2 2 1 ) can only be made in one way, so that 
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l^- 2 2 i ) + Ce L . 
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To find the coefficient C we just count terms: NReL-i 
has (JV + 1)!/(JV - L + l)l(N - 1)! terms, \A) has 
Nl/(N-L + l)l(L-2)l terms, and e L has N\/(N-L)\L\ 
terms. This leads to the result C — L, and hence 



\A) = NReL-i - Le L . 
Combining this with Eq. (Elf) gives 



(24) 



He L = 



N' A - I 1 + - I JV 



N(N — L + l) _ 
e L + — z L Re L -i- 



(25) 



(5) Operating H onto Sm 

If we now operate H onto e-M using Eq. jl7| ) to write cm 
in terms of ex, and R, and using Eq. ( p5| ) to operate H 
onto ez, we get after a lot of tedious algebra, 
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N 2 



M 
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(26) 



Note that we have to compare the coefficients of three 
types of term separately: (i) e^, (ii) R M ^ L e^ for 2 < 
L < M, and (hi) R . The proof can be simplified a little 
if we use the subspace property. Introduce a projection 
operator, P, which removes any term containing a factor 
of R. Operating P onto Eq. p7|) gives Pcl = e.L, whilst 
operating P onto Eq. (Ea) gives Eq. (Efl). 
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IV. DISCUSSION AND CONCLUSIONS 

We have considered N bosons in a 2D harmonic po- 
tential interacting via repulsive delta-function potentials 
and with fixed total angular momentum L < N. Within 
the "lowest Landau level" approximation, we have ana- 
lytically shown that the elementary symmetric polyno- 
mial of coordinates relative to the center of mass, &l is 
an exact eigenstate of this Hamiltonian with eigenvalue 
N(N — L/2 — 1). Extensive numerical analysis shows that 
this state is actually the ground state. This is not sur- 
prising since in the special case L = N this is just the 
one-vortex state discussed in Ref. || which is expected 
to be the ground state by analogy to superfluid 4 He. One 
can also see that there is a sense in which distributes 
the angular momentum equally between particles subject 
to the subspace property. We have attempted to prove 
analytically that is the ground state, but have so far 
failed. The situation is much harder than in the proof we 
presented in Ref. || to show that the eigenstate with the 
largest eigenvalue corresponds to all angular momentum 
being in the center of mass motion. The main problem 
is that the eigenvector for the smallest eigenvalue has 
components of both signs to reduce its energy, and frus- 
tration results between any set of three basis states that 
are connected by H: the difficulty is essentially that of 
a quantum antiferromagnet compared to a ferromagnet. 
For completeness we hope that the state considered in 
this paper will be analytically shown to be the ground 
state in the future, but from numerics there can be little 
doubt that it is. 
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